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Abstract 



>. 

, B s — ► p(uS)K* are useful to determine the B s distribution amplitude, as well as constrain the 

. CKM phase angle a. We study these decays within the Perturbative QCD approach(PQCD). 



In this approach, we calculate factorizable, non-factorizable, as well as annihilation diagrams. 



We find the branching ratio for B s — ► p + K* is at the order of 10 5 , and there's large direct 



CP violation in B S {B S ) -> p°K*°(uK* ). Our predictions are consistent with those from other 



methods and current experiments. 



1 Introduction 



X 

Exclusive nonleptonic B decays have provided a fertile field to investigate the CP violation 



and search for new physics. The hadronic matrix elements of the effective operators play a key 
role in the study of B meson decays, but it is difficult to calculate them precisely due to the long 
distance QCD dynamics. The factorization approach (FA) p] [2] based on the color transparency 
mechanism has been applied to many decay modes, and it works well in many channels. But it 
suffers from some problems such as infrared-cutoff and scale dependence. To solve these problems 

and make more accurate predictions, the perturbative QCD approach (PQCD) [H1IU El, the QCD 
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improved factorization (QCDF) [3IH] as well as the Soft-collinear effective theory (SCET) 9 have 
been developed in the recent years. 

PQCD is based on kj< factorization theorem |10UllllT2] . The decay amplitude is factorized into 
the convolution of the mesons' light-cone wave functions (see Appendix A), the hard scattering 
kernels and the Wilson coefficients, which stand for the soft, hard and harder dynamics respectively. 
The transverse momentum is introduced so that the endpoint singularity which will break the 
collinear factorization is regulated and the large double logarithm term appears after the integration 
on the transverse momentum, which is then resummed into the Sudakov form factor. The formalism 
can be written as: 



M ~ J dxidx2dx 3 b 1 dbib2db2b 3 db3Tr[C(t)<$>B(xi,bi)$ K «(x2,b 2 )®p{x3,b 3 ) 

Hixi^mi^e- 3 ^}, (1) 

where the bi is the conjugate space coordinate of the transverse momentum, it denotes the transverse 
interval of the meson, t is the energy scale in hard function H. The jet function St(xi) comes from 
the summation of the double logarithms In 2 Xi near the endpoint, called threshold resummation 
|10l I13j . The factorization theorem guarantees the infrared safety and the gauge invariance of the 
hard kernel and has been proved to all order of a s |14j . 

Many hadronic two body B decays have been studied in PQCD approach [Hj l6| IT5 | I16j. Most 
predictions are consistent with the current experiments. The B s decays are important to extract 
CKM phase angles and study the CP violation. As B s meson is not in the energy scale of the high 
luminosity B factories SLAC and KEK, it is more difficult to be produced and measured now. We 
can study the B s decays more precisely in the very near future with the increase of luminosity at 
TEVATRON and the upcoming Large Hadron Collider (LHC). 

B s meson is different from B meson due to the heavier strange quark (compare to u, d quark) 
which induces the SU(3) symmetry-breaking effect. This effect is considered to be small and the 
distribution amplitude of B s meson(given in the following formula) should be similar to that of the 
B meson, 

l,xM Bs . 2 uj 2 Bs b 2 



b Bs {x) = N Ba x 2 (l- x) 2 exp 



The upper limit of the B s — > ttK branching ratio is 7.5 x 10 -6 |17j . which constrain the parameter 
ujb s to a lower limit of about 0.5 ^Sj. Moreover, in order to fit the branching ratio measured 
in the B s — > cfxf) decay |19j . we constrain ub s to about 0.55 [201; then we can see that the SU(3) 



symmetry-breaking is not negligible. Here we integrate out the variable b and show the distribution 
amplitude of B and B s meson in Fig.l. We can see that the peak point of the curve of B s meson's 
distribution amplitude prefers a larger x (x denotes the momentum fraction carried by the light 
quark) region comparing to the B meson. This is consistent with the fact that the s quark much 
heavier than the d (u) quark, should carry more momentum. Later in this paper, we will see that 
the branching ratios of B s decays are very sensitive to this parameter. If measured by experiments, 
radiative leptonic decays of B s meson can provide information of this parameter |2J . 



1.3 r- 
1.2 - 




0.50 



Figure 1: distribution amplitude as a function of light quark momentum fraction x (The wave function comes 
to zero when x £ {0.5, 1.0} so we don't show it here) 

In this paper, we study B s — > p(uj)K* decays in the PQCD approach. Hopefully the branching 
ratio is not too small and can be detected by the TEVATRON or LHCb experiments, then it may 
allow us to determine the B s distribution amplitude and SU(3) breaking effects with much more 
precision. Moreover, we can also constrain a with fewer pollution from this channel. 

2 Calculation and Numerical analysis 

We use the effective Hamiltonian for the process B s — > p{u)K* given by |22j 

n eff = ^§ |k [CxWom + c 2 ( m )o 2 (m)] - v t ^q^o^V)} , (3) 
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where V u = V^Vub, Vt = V^tfo, Ci(m) are the Wilson coefficients, and the operators are 

0\ = (diqj) V -A(qjbi)v-A, 2 = (diqi) V -A{qjbj)v-A, 

O3 = (dibi)v-A^2(q~jqj)v-A, C>4 = (dibj)v-A^2(<ijqi)v-A, 
q q 

O5 = (dibi)v-A^2(q~jqj)v+A, Oq = (dibj) V -A^2(qjqi)v-A, 

q 1 
2 _ 2 _ 

°7 = -(dibi) V -A^2e q (qjqj) v+ A, O s = -(dibj) V -A^2e q (qjqi)v+A, 

<? 1 
2 _ 2 _ 

Og = -{dibi) V -A ^2 e <?(<7i<?i)v-A, O10 = ^{dibj) V -A ^ e^q^v-A- (4) 

Here i and j stand for SU(3) color indices. 
The decay width for these channels is : 

g^IpI 



16ttM 2 



M a] M u (5) 

L B a=L ,T 

where p is the 3-momentum of the final state mesons, |p| = ^^(l — r 2 K * — r "^^)i an d ^x*(p,w) = 
m K*(p,uj)/ m B s - M a is the decay amplitude , which will be calculated later in PQCD approach. The 
subscript a denotes the helicity states of the two vector mesons with the longitudinal (transverse) 
components L(T). According to Lorentz structure analysis, the amplitude can be decomposed into: 

M" = Ml Ml + M 2 B M N el{a = T) ■ e* 3 (a = T) + iM T e^ p(T e%* efPpf . (6) 

We can define the longitudinal Hq, transverse H± helicity amplitudes 



H = M 2 B M L , H± = Ml s M N =F M\, vV 2 - 1M T , (7) 
where r' = ,/ 2 'w — . They satisfy the relation 

M aj M a = \H \ 2 + \H + \ 2 + \H_\ 2 . (8) 

a=L,R 

The leading order diagrams in PQCD approach are shown in Fig. 2. The amplitudes for B s — > 
pK*° and B s — ► pK*°&re written as 

m h = v:t h -v;p h , (9) 

M H = V u T H -V t P H , (10) 

respectively, where the subscript H denote different helicity amplitudes L, N and T, and Tjj and 
Ph are the amplitudes from tree and penguin diagrams respectively. The detailed formulae of Tjj 
and Ph are similar to those in B — > K*K* [2S] and B — ► <pK* JH], so we will not show them here. 




Figure 2: Leading order Feynman diagrams for B s — > p(u)K* 

The parameters used in our calculations are: the Fermi coupling constant Gf = 1.16639 x 
W- 5 GeV- 2 , the meson masses M Bs = 5.37GeV, M K * = 0.89GeF, M p{w) = 0.77GeV 24 , the 
decay constants f K * = 0.217GeV, f%, = 0.16GeV, f p = 0.205GeV, f p = 0.155GeF, / w = 
0.195Gey, /J = O.UGeV (23], the central value of the CKM matrix elements a = 95°, \V td \ = 
0.0075, \Vtb\ = 0.9992, \V ud \ = 0.9745 (21, \V ub \ = 0.0047 : 26 and the B s meson lifetime 
t Bs = 1.461ps (2H- 

As we mentioned before, the decay i? s — > p{uj)K* can be used to determine the B s meson 
wave function parameter wg s , or u>b 3 can influence our predictions of B s — > p(u>)K* decay. So we 
show the results in Table 1 according to 3 different values of wb s . From the table we can easily 
find out the averaged branching ratio for B S (B S ) — ► p^K*^ is much larger than the other two, for 
B S (B S ) — > /r^iT*^ involve large Wilson coefficient C(2) + C(l)/3 for the factorizable part while 
the other two (B S (B S ) -> p°K*°(K* ) and B S (B S ) wif *°(K* )) involve a much smaller Wilson 
coefficient C(l) + C(2)/3 (color-suppressed) for the factorizable part of the emission diagram. As a 
result the first one is tree dominated and has a large branching ratio and small direct CP asymmetry. 
While referring to the other two, the contributions from penguin and tree diagrams are at the same 
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Table 1: Branching ratio, polarization fraction and direct CP asymmetry dependence on lob s 



Channel 




BR(10" 6 ) 


\A \ 2 


|^l|| 2 


l^l 2 


Adir 




0.50 


0.41 


0.38 


0.28 


0.34 


-0.63 


B S (B S ) -» p°K*°(K* ) 


0.55 


0.34 


0.41 


0.27 


0.32 


-0.67 




0.60 


0.30 


0.44 


0.26 


0.30 


-0.70 




0.50 


0.56 


0.33 


0.31 


0.36 


0.56 


B S (B S ) -> luK*°(K*°) 


0.55 


0.47 


0.35 


0.30 


0.35 


0.58 




0.60 


0.40 


0.37 


0.29 


0.34 


0.60 




0.50 


16 


0.92 


0.04 


0.04 


0.10 


B s (B s )^p ± K*^ 


0.55 


12 


0.92 


0.04 


0.04 


0.12 




0.60 


10 


0.92 


0.04 


0.04 


0.13 



order {Zh — 0.5 ~ 1-5), hence we can expect a large direct CP asymmetry from eas.(|13|). 

The polarization fraction difference of these channels are also due to that the main contribu- 
tion of each channel comes from different topology. B S (B S ) — > p^K*^ is tree dominated. The 
main contribution comes from the factorizable part of the emission diagram, where transverse 
polarization amplitude is suppressed by a factor r? ~ (0.77/5.37) 2 (see formulas in 23j), so the 
longitudinal polarization dominates and contributes more than 90% of the total branching ratio. 
But in B S (B S ) -> p°K* (K*°) and B a (B a ) ujK*°(K*°) decays, tree emission (factorizable) dia- 
gram contribution is suppressed due to the cancellation of Wilson coefficients C\ + C2/3. The left 
dominant contribution is the non-factorizable diagrams of tree operators and penguin diagrams, 
both of these contributions equally contribute to longitudinal and transverse polarizations. The 
transverse polarization is not suppressed in those cases, therefore numerically we get a small longi- 
tudinal fraction of about 0.4. This similar situation is also found in B — > pp(uj) decays |27| . which 
are related by SU(3) symmetry to our B s — > pK* decays. 

To extract the CP violation parameters and dependence on CKM phase angle a of these decays, 
we rewrite the helicity amplitudes in 1)9110(1 as the functions of a: 

m+ = v:t h -v;p h 

= V*T H (1 + Z H e^ a+5 ^) (11) 
M H = V u T H -V t P H 

= V u T H (l + Z H e i (~ a+SH '>) (12) 
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Figure 3: Average branching ratios when ub s was set from 0.50 to 0.60: (a)B s (B s ) — > p°K*°(K*°); (b) B 3 (B S ) — » 




0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 



Figure 4: Average branching ratios for B S (B S ) — > p ± K*^ 1 as a function of a when o;.b s was set from 0.50 to 0.60 

where Zh = \ V* /V*\\Ph /Th\, and 5 is the relative strong phase between tree (T) and penguin (P) 
diagrams. Here in PQCD approach, the strong phase comes from the nonfactorizable diagrams and 
annihilation diagrams. This is different from Beneke-Buchalla-Neubert-Sachrajda 7 approach. In 
that approach, annihilation diagrams are not taken into account, strong phases mainly come from 
the so-called Bander-Silverman-Soni mechanism 28 . As shown in [Sj, these effects are in fact 
next-to-leading-order (a s suppressed) elements and can be neglected in PQCD approach. We give 
the averaged branching ratios of B S (B S ) — > p° K*° (K*°) and B S (B S ) — > u>K*°(K*°) as a function 
of a in Fig. 3, and the averaged branching ratios of B°(B°) — > p^K*^ in Fig. 4. 



Using Eas. 1)111121) . the direct CP violating parameter is 

Ad „ = \M\ 2 -\M\ 2 
CP \M\ 2 + \M\ 2 

—2sina{T 2 sinSL + 2T^sinSN + 2T^sinST) 
T|(l + Z\ + 2Z L cosacos5 L ) + 2 J2i =N ,T T i { 1 + z ? + 2Z i cosacos5 i ) ' 

Notice the CP asymmetry for these channels are sensitive to CKM angle a, we show the direct CP 
asymmetry as a function of a in Fig. 5. It is easy to see that the B s — > p°K*° and loK*° have large 
direct CP asymmetries up to 50%, with a relative minus sign. On the other hand, the B s — > p^K * =F 
decay has small direct CP asymmetry due to only one large tree contribution in this decay. The 
uncertainty shown at this table is only from the B s meson wave function parameter dependence. In 
fact, since CP asymmetry is sensitive to many parameters, the line should be more broadened by 
uncertainties. The mixing induced CP asymmetry is complicated and requires angular distribution 
study, similar study may be found in |29j . 
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Figure 5: Direct CP asymmetry A^ r p as a function of a: (a) B S (B S ) — > p°K*°(K*°) (lower line); 
(b) B S (B S ) -> loK*°(K* ) (upper line); (c) B S (B S ) -» p ± K* T (middle line). 



At last, if we compare our predictions with those of naive factorization [HO] 

BR{B S (B S ) -» p°K*°(K* )) = 5.5 x 1(T 7 , (14) 

BR{B S {B S ) -> uK*°(K* )) = 6.0 x 10" 7 , (15) 

BR(B S (B S ) -> ^K*^) = 1.7 x 10~ 5 , (16) 
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and ones of QCDF [3T] 



BR(B S (B S ) -> p { 



5.3 x 10 



7 



(17) 



BR{B S {B S ) 



uK*°(K* )) 



3.1 x 10 



7 



(18) 



BR{B S {B S ) p ± K* T ) 



1.8 x 10 



-5 



(19) 



We can see that they are consistent. It should be noticed that the branching ratios in FA and 
QCDF strongly depend on form factors. While in PQCD, the branching ratios and form factors 
depend on wave functions, especially the B s meson wave function. Nowadays, very few B s meson 
decays have been measured, so we can only give rough constraints on the parameters from other 
channel and permit large errors. More experimental data can help to constrain the form factors and 
wave functions, then we can give more precise predictions and the different methods can be tested 
by the experiments. Although similar results are got by different methods for branching ratios, 
the polarization fractions are quite different. The QCDF and naive factorization give only several 
percent transverse polarization for all three decay modes [3^, while Table 1 shows large transverse 
contribution for B s — > p°(uj)K* decays in our PQCD approach. The direct CP asymmetry are not 
given in ref.jSJ^ but they probably also differ from PQCD approach as it happened in B — > tttt 
and Kit case |32| . 

The numerical results shown here are only leading order ones. For the tree dominated channel 
B s — > p K*^ , the leading order diagrams should give the main contribution. But for the other two 
decays, with a branching ratio as small as 10 -7 , the next-to-leading order and power suppressed 
contributions should not be negligible, the results may suffer from large corrections when the next 
to leading order corrections are included |33| . 

Current experiments [2l] only give the upper limit for the decay 



More data are needed to test our calculations. 

3 Summary 

In this paper we calculate the branching ratios, polarization fraction and CP asymmetries of 
B s — > p(uj)K* modes using PQCD theorem in SM. We perform all leading order diagrams to next 
to leading twist wave functions. We also study the dependence of their averaged branching ratios 



BR(B S -> p°K*°) < 7.6 x 10 



-4 



(20) 
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and the CP asymmetry on the CKM angle a. At last we compare our predictions with values from 
other approaches. 
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A wave function 

For longitudinal polarized K* meson, the wave function is written as 



1 



(21) 



and the wave function for transverse polarized K* meson reads 



1 



e» 2T P 2 p n<Lr K *(x)]. 



(22) 



1€ liu polSl' 
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The K* meson distribution amplitudes up to twist-3 are given by ref.[3I] with QCD sum rules. 

4> K *{x) = ^MLx(l-x)[l + 0.57(1 -2x) + 0.07C 2 /2 (l-2x)], (23) 
V 2iV c 

<^,(x) = {0.3(1 - 2x)[3(l - 2x) 2 + 10(1 - 2s) - 1] + 1.68C 4 1/2 (1 - 2x) 

+0.06(1 - 2x) 2 [5(l - 2x) 2 - 3] + 0.36 {1 - 2(1 - 2a:) [1 + ln(l - x)]}} , (24) 

fT 

<f) s K * (x) = ^= {3(1 - 2x) [l + 0..2(1 - 2x) + 0.6(10x 2 - lOx + 1)] 

-0.12x(l - x) + 0.36[1 - 6x - 2 In (1 - x)]} , (25) 
f K . (x) = ^=x(l - x)[l + 0.6(1 - 2x) + 0.04C 2 3/2 (1 - 2x)], (26) 

(x) = -%= { ~ [1 + (1 - 2x) 2 + 0.44(1 - 2x) 3 ] 



2^2lVc l4 L 

+0.4C 2 /2 (1 - 2x) + 0.88C 4 1/2 (1 - 2x) + 0.48[2x + ln(l - x)]} , (27) 

fT 

4>K*{x) = 7^— {3(1 - 2s) [1 + 0.19(1 - 2x) + 0.81(10x 2 - lOx + 1)] 

-1.14x(l - x) + 0.48[1 - 6x - 2 ln(l - x)]} , (28) 
where the Gegenbauer polynomials are 

cho = \(ze -i), (29) 

cl(O = ^(35C 4 -30^ 2 + 3), (30) 

C 2 f (0 = §(5e 2 -1). (31) 

For p and u> meson, we employ p° = -^(uu — dd) and u> = -^(uu + dd). Their Lorentz structures 
are similar to K* meson, the distribution amplitudes are the same for p and uj and given as |34j : 

<j> p {x) = _^x(l-x)[l + 0.18C 2 3/2 (l-2x)] , (32) 

#(*) = ^-^{3(l-2x) 2 + 0.3(l-2x) 2 [5(l-2x) 2 -3] 

+0.21[3 - 30(1 - 2x) 2 + 35(1 - 2x) 4 ]} , (33) 
3f T 

<j> s p (x) = _^=(l-2x)[l + 0.76(10x 2 -10x + l)] , (34) 

4> T p (x) = -^x(l-x)[l + 0.2C 2 3/2 (l-2x)] , (35) 



<f>;(x) = -^|r [1 + (1 _ 2x) ^ +0 .24[3(l-2x) 2 -l] 



+0.12[3 - 30(1 - 2x) 2 + 35(1 - 2x) 4 ] j , (36) 
r p (x) = ^^(l-2x)[l + 0.93(10x 2 -10x + l)] . (37) 
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